Finite-distance singularities in the tearing of thin sheets 
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We investigate the interaction between two cracks propagating in a thin sheet. Two different 
experimental geometries allow us to tear sheets by imposing an out-of-plane shear loading. We find 
that two tears converge along self-similar paths and annihilate each other. These finite-distance 
singularities display geometry-dependent similarity exponents, which we retrieve using scaling ar- 
guments based on a balance between the stretching and the bending of the sheet close to the tips 
of the cracks. 
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Thin sheets are widespread in Nature and technology. 
Examples include insect wings, leaves or tectonic plates, 
and graphene, conducting layers or metallic roofs, respec- 
tively. Strikingly, these objects feature two types of en- 
ergy focusing, at a sharp fold as in crumpled paper [l| or 
at the tip of a crack as in a torn sheet of paper. More 
generally, understanding and predicting the propagation 
of a crack in a brittle material yield a central challenge 
in fracture mechanics [2, 3]. A body of research has been 
performed on the path of a crack submitted to in-plane 
tensile (mode I) or shear (mode II) loading, see (j-[6| for 
recent references. An emerging consensus is that a crack 
generally propagates along a path such that mode II van- 
ishes and the material opens in a tensile mode (principle 
of local symmetry). Few studies considered out-of-plane 
shear (mode III) loading, see 0, 0. Even fewer inves- 
tigations addressed mode III crack propagation in thin 
sheets, focussing on ductile materials [9l4llj|, on sheets 
adhering to a substrate [HI, [l3j , or on cutting with blunt 
indentors pH EH, [H|. This class of experiments point 
out the coupling between the bending of the sheet and 
in-plane stretching leading to the opening of cracks. 

In a preliminary report [16|, we studied the propaga- 
tion of one or two cracks in a brittle thin sheet, using the 
so-called trousers test configuration [9] . The results that 
we obtained indicate that, although a large scale mode 
III is imposed, the material locally breaks in the tensile 
mode I. As a consequence, the interfaces opened by the 
crack are not perpendicular to the sheet surface and are 
inclined by 45 degrees instead. Here we also consider a 
peel-like configuration and we focus on two tears propa- 
gating quasi-statically. We find that they converge along 
self-similar paths, with characteristic exponents corre- 
sponding to each type of loading. The topological change 
that occurs when they annihilate each other is reminis- 
cent of finite-time singularities in the breakup of liquid 
droplets and jets [13, [H| , which suggests the terminology 
of finite-distance singularity. Our power-law crack paths 
differ from the exponential shapes predicted in [l9j; we 
account for the exponents that we observe using scaling 
arguments inspired by energy localization along ridges in 
crumpled sheets [20|. Therefore our experiments appear 
as a non-trivial combination of the two types of energy 
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FIG. 1: The trousers configuration: side (a) and top (b) 
schematics of the experiment. The sheet is pulled from the 
three strips as shown by the arrows. The strips are deflected 
by freely rotating cylinders and then rolled over two cylinders 
entrained by a continuous motor, (c) After the tears have con- 
verged, the central strip has a tongue-like shape, described by 
a curve y(x). 



focusing known in thin sheets. 

In experiments, we used thin films of bidirectional 
polypropylene of thickness h = 30, 50 and 90/im. The 
Young's modulus and fracture toughness were measured 
as 2.2±0.4 GPa and 2.6±0.3MPa.m 1 / 2 , respectively. The 
films were found to have isotropic mechanical properties 
within these uncertainties. All experiments were per- 
formed quasi-statically, with crack velocities in the range 
0.05-1.5mm.s _1 . At ambient temperature and in this ve- 
locity range, the fracture process is brittle for polypropy- 
lene. The crack paths were digitized using a scanner, with 
a resolution of about 10/im. 

The first setup is inspired by the trousers test used to 
characterize ductile sheets [9[ and mimics the common 
way one tears a sheet of paper in three pieces. It was 
described in [l6| ; we recall here the important features of 
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the setup. In a very long sheet of width W, two cracks are 
initiated by two cuts positioned symmetrically with re- 
spect to the central axis of the sample, so that three strips 
are created at one end of the sheet (Fig. [1]). The sheet 
is pulled from the three strips using cylinders for the en- 
trainment. The main control parameters are the width of 
the sample W, the initial distance between the two cracks 
w, and the distance d between the freely rotating cylin- 
ders that deflect the strips. We found that the results 
presented here are independent of d. When the tears are 
propagated the distance between their tips, Y, decreases 
from w to 0, as they converge and eventually annihi- 
late each other. The topology has changed: the sheet 
is split in two parts, and the central strip has detached 
into a tongue-like shape, as shown in Fig. QJ. When the 
ratio between the initial width of the central strip and 
the width of the sample is smaller enough (w < 0.1 W), 
experiments are reproducible and the results are found 
to be independent of W. If the central strip is larger 
(w > 0.1 W), experiments are less reproducible and the 
shape of the tear is not always symmetric. This might 
be ascribed to the sensitivity of propagation to imperfec- 
tions in the parallelism of the cylinders' axes. 

We analyzed the crack paths by digitizing the tongue- 
like central strip. The width y of the strip is defined in 
Fig. and shown in Fig. [2^ as a function of the dis- 
tance x to the point where the two tears converged. An 
important observation is that all realizations are super- 
imposed (for various values of w, W, and d), indicating 
that the geometry of the setup is unimportant for path 
selection. We found that, in the range 5 < x < 300mm, 
the curves y(x) are well-described by a power-law of expo- 
nent 0.64 ±0.06 and a prefactor (0.77 ±0.2 and 1.2 ±0.2) 
that is larger for the largest value of thickness. Below 
x = 5mm, y(x) seems close to being linear. We then 
sought the origin of this transition and made the follow- 
ing observations. After tearing, far from the tip, the strip 
recovers its flatness indicating a brittle fracture process, 
while it is permanently curved in a small region close 
to the tip, which is a signature of plastic deformations. 
Post-mortem analysis of the fracture surface shows that 
it is generally inclined by an angle /? = 45° with respect 
to the surface of the film, in agreement with our previous 
observations on a single crack [l6[. However, a transi- 
tion occurs for x ~ 5mm: this angle decreases down to 
the value /3 = 34° at the merging point. This transition 
seems to be related to the appearance of plastic deforma- 
tions. If (3 decreased to zero, the fracture surfaces would 
be parallel to the sheet surfaces, and the configuration 
would resemble the peeling of an adhesive strip, in which 
two notches lead to triangular strips [12|. This remark 
might account for a linear behavior of y(x) for small x. 
Finally, we non-dimensionalized the profile y(x) by us- 
ing the thickness h as a unit (Fig. [2fc), which yields a 
fair collapse of the data. However, we were unable to 
obtain reproducible results with the smallest thickness, 
h = 30/im, for which the crack paths are very sensitive 
to the forcing. As a consequence, the range in thickness 
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FIG. 2: Crack paths in the trousers configuration, (a) Shape 
of the central strip: width y as a function of longitudinal coor- 
dinate x, see Fig.[Tt for definitions. The widths for a thickness 
of 90/xm were multiplied by 100 for clarity. The colors of the 
symbols correspond to different realizations. Power laws of 
exponents 2/3 and 1 are shown for comparison, (b) Angle f3 
of inclination of the fracture surface with respect to the sheet 
surface (defined in inset) as a function of the longitudinal dis- 
tance x, suggesting a transition around x = 5mm from 45° to 
lower values. The initial distance between the two tears was 
w — 20mm and the thickness had the value h — 50/mi. (c) 
The curves in (a) were made nondimensional using the sheet 
thickness h as a unit. Black and red symbols correspond to 
h = 50/xm and 90/im, respectively 
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FIG. 3: The peel-like configuration, (a) Schematic of the 
experiment. Two notches are initiated in a rectangular long 
sheet clamped along its lateral boundaries. The central strip 
is pulled horizontally at constant velocity until it detaches 
when the two tears converge, (b) Tongue-like shape of the 
detached strip, defined by y(x). 



is too narrow to rule out other small-scale characteristic 
lengths. 

In a second step, we considered a less symmetric, 
peel-like experimental configuration, which was inspired 
from [12]. A long sheet of width W is clamped along 
its lateral boundaries to a thick wooden plate, using nar- 
row adhesive tapes (Fig. [3^). Two parrallel notches are 
initially made at a distance w one from an another at 
one end of the sheet. The central strip is pulled hori- 
zontally, so that the two tears propagate quasi-statically. 
In contrast with the first configuration, the distance be- 
tween the pulling point and the crack tips increases, but 
this macroscopic length appears to be unimportant in 
the following results. As the tears advance, the distance 
between the two tips decreases from w to 0, when they 
annihilate each other and the central strip detaches. The 
resulting shape y(x) (Fig. [3b) is qualitatively similar to 
those of the first setup. The shapes of various detached 
strips is shown in Fig. H^i. Again, all realizations are su- 
perimposed for a given value of thickness, indicating that 
the geometry of the setup is unimportant for path selec- 
tion. Over two orders of magnitude, the curves y(x) are 
well described by a power-law 0.77 ±0.05 and a prefactor 
(0.82 ±0.06, 1.2 ±0.2, and 1.53 ±0.2) that increases with 
thickness. As in the first setup, the tip of the central 
strip undergoes plastic flow, but over a smaller length 
(~ 1mm). However, we have not found any signature 
on the crack paths. Finally, we non-dimensionalized the 
profile y(x) using the thickness h as a unit (Fig. HJd), 
which provides a satisfying collapse of the data. 

Overall, the propagation of the two tears leads to a 
topological change such that the cracks annihilate each 
other and the central strip detaches. The crack paths 
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FIG. 4: Crack paths in the peel-like configuration, (a) Shape 
of the central strip: width y as a function of longitudinal coor- 
dinate x, see Fig. [3)3 for definitions. The widths for the values 
of the thickness 30/rni and 90/rni were divided and multiplied 
by 10, respectively, for clarity. The colors of the symbols cor- 
respond to different realizations. A power laws of exponents 
3/4 is shown for comparison, (b) The curves in (a) were made 
nondimensional using the sheet thickness h as a unit. Green, 
black and red symbols correspond to h = 30/xm, 50/xm, and 
90/im, respectively 



seem to follow well-defined power-laws. In order to ex- 
plain this behavior, we consider the ridge joining the two 
tips, which is one form of energy localization in thin 
sheets [1] as the sheet is highly bent. Guided by the 
principle of maximal release of energy for the selection 
of a crack path, we postulate that the two cracks follow 
the stress field induced by the presence of this ridge and 
we investigate how this stress field decays away from the 
ridge. This ridge of length Y imposes a normal displace- 
ment a that is felt over a distance X in the region ahead 
of the cracks (see Figs. QJi,b and |3ji). Following [20|, we 
estimate the length X by balancing the stretching en- 
ergy S s and the bending energy E\> in a region of area 
S ~ XY. Assuming Y <C X, the strain s is distributed 
over the larger dimension of this region so that it can 
be estimated as s ~ (a/X) 2 , yielding a stretching energy 
S s ~ Ehs 2 S. The main curvature c is along the shortest 
direction, so that c ~ a/Y 2 , yielding Sh ~ Eh 3 c 2 S. As a 
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consequence, £ s ~ £^ corresponds to 

Y~x{h/ a y'\ (i) 

which is consistent with the assumption that Y <C X as 
long as a > L 

In the case of the trousers configuration, the normal 
displacement a corresponds to a slope 6 felt over a dis- 
tance X/2, so that a ~ X6. In experiments, we observed 
this slope ~ 0.1 to be roughly constant, but we do not 
have a theoretical argument for its selection. Replacing 
in equation (pQ), we obtain 

Y ~ (h/Of^X 2 ^. (2) 

Due to the weak dependance on we expect prefactor 
of order 1 for Y/h ~ (X/h) 1 / 3 . Following the postulate 
that the cracks follow the shape of the region where the 
ridge is felt, this scaling should predict the paths. Indeed, 
it yields a good description of experimental results in 
Fig. [2t, except for the linear behavior for small x. 

In the peel-like configuration, the normal displacement 
is given by the width of the ridge boundary layer [20| . 
a ~ /i 1 / 3 X 2 / 3 . A substitution in equation (pp) leads to 

Y ~ h^X 3 ^. (3) 

This is also in good agreement with the experimental 
crack paths shown in Fig. 0)3. 

Thus, we are able to account for the power-law behav- 
ior in experiments by estimating the decay length of the 



stress field penetration in the sheet. This behavior is in- 
dependent of macroscopic parameters but the exponents 
depend on the symmetries of the experimental configu- 
ration, as in the finite-time singularities occurring in the 
breakup of liquid droplets and jets [171. 118] . This similar- 
ity, such that the coordinate along the axis of the sheet 
replaces time, prompted us to use finite-distance singu- 
larities to describe our findings. Our results differ from 
the exponential shapes predicted in [19], possibly because 
our sheets are an order of magnitude thinner than in their 
simulations. Indeed, in the case of the propagation of a 
single crack, the numerical results of Cohen et al. [19| 
get closer to our experiments when they decrease the 
thickness. Future theoretical work should address this 
discrepancy, as well as other features that we have not 
accounted for, such as the selection of the slope 9 or the 
transition to a three-dimensional fracture. Although the 
coupling between bending and stretching seem to imprint 
the stress field and guide the tears, the principles under- 
lying propagation of cracks in out-of-plane shear loading 
are still to be established. 
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